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Abstract 

a: 

We discuss models of neutrino masses that lead to a large mixing angle for atmospheric neutrino 
oscillations. In particular we study a mechanism where a simple texture for the Dirac matrix 
leads to the observed pattern of mixings for a continuous range of Majorana matrices in a 
reasonably natural way Both nearly maximal and small mixing solutions for solar neutrinos 
are compatible with this model. Possible dynamical realizations are discussed and a detailed 
example in terms of horizontal U(l) charges is presented. 
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Following the experimental results from Superkamiokande JIJ a lot of attention has been 
devoted to the problem of a natural explanation of the observed nearly maximal mixing angle for 
atmospheric neutrino oscillations. It is possible that also solar neutrino oscillations, if explained 
by vacuum oscillations, occur with a large mixing angle |2| . Large mixing angles are somewhat 
unexpected because the observed quark mixings are small and the quark, charged lepton and 
neutrino mass matrices are to some extent related in Grand Unified Theories. In a previous 
paper || we have given a first discussion of this problem. Within the framework of the see-saw 
mechanism we focussed on the interplay between the neutrino Dirac and Majorana matrices 
which is necessary in order to generate maximal mixing Q . In the present note we continue the 
study of this issue. We start by giving a general formulation of the problem and by specifying 
our working assumptions. Among our previously proposed strategies for a natural explanation 
of maximal mixing, we select and further elaborate on one which we find particularly plausible 
and appealing. For our mechanism to work we need that, in analogy with what is observed 
for quark and leptons, one mass eigenvalue is dominant for neutrinos (for both the Dirac and 
the effective light neutrino mass matrices). Also, in zeroth approximation, we need that the 
neutrino Dirac mass matrix, in the basis where charged leptons are diagonal, has only two large 
entries of comparable magnitude. We then show that either single or double nearly maximal 
mixing occur automatically for whatever Majorana mass matrix chosen in a large domain of 
parameter space, without need of fine tuning. We discuss models where these conditions can be 
realized and give numerical examples that fit the present data. We finally compare ours with 
other possible mechanisms. 

We start by assuming only three flavours of neutrinos that receive masses from the see-saw 
mechanism |5[] and allow all possible hierarchical patterns for the neutrino mass eigenvalues m^. 
For reasons of simplicity, we consider the simplest version of the see-saw mechanism with one 
Dirac, mp, and one Majorana, M, mass matrix, related to the neutrino mass matrix m v , in 
the basis where the charged lepton mass matrix is diagonal, by 



m, 



mF D M~ x m D . (1) 



As well known this is not the most general see-saw mechanism because we are not including 
the left-left Majorana mass block. Maximal atmospheric neutrino mixing and the requirement 
that the electron neutrino does not participate in the atmospheric oscillations, as indicated by 
the Superkamiokande [|1| and Chooz data, lead directly to the following structure of the U /j 
(f=e,/i,r, i=l,2,3) real orthogonal mixing matrix, apart from sign convention redefinitions (here 
we are not interested in CP violation effects: all matrices are taken real) 



U 



fi 



C S 

s/y/2 c/y/2 -1/V2 
s/y/2 c/y/2 +l/y/2 



(2) 



This result is obtained by a simple generalization of the analysis of ref. Q (also discussed in 
ref. ||) to the case of arbitrary solar mixing angle (s = sin 9 sun , c = cos 9 sun ): c = s = 1/ y/2 for 
maximal solar mixing (e.g. for vacuum oscillations sin 2 2# sun ~ 0.75) , while sin 2 29 sun ~ 4s 2 ~ 
5.5 • 10~ 3 for the small angle MSW || solution. The vanishing of U e3 guarantees that v e does 
not participate in the atmospheric oscillations and the relation \U^\ = \U T s\ = I /y/2 implies 
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maximal mixing for atmospheric neutrinos. Note that we are assuming only two frequencies, 
given by 

A SMn (xm 2 2 - m 2 , A atrn ocnij- m 2 2 (3) 

The effective light neutrino mass matrix is given by m v = Um,nagU T with m^ag = Diag\m\, m 2 , m 3 ] 
For generic s one finds 

^2e 5 5 
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with 



e = (mic 2 + m 2 s 2 )/2 , 5 = (mi — m 2 )cs/v2 



£2 



(mis 2 + m 2 c 2 )/2 



(4) 



(5) 



With respect to the notation of ref. || we have reinstated the normalization factor m 3 /2 in 
m v . We see that the existence of one maximal mixing and U e3 = are the most important 
input that leads to the matrix form in eq. (§,^J). The value of the solar neutrino mixing angle 
can be left free. While the simple parametrization of the matrix U in eq. (|2p is quite useful to 
guide the search for a realistic pattern of neutrino mass matrices, it should not be taken too 
literally. In particular the data do not exclude a non- vanishing U e ^ element. In most of the 
Superkamiokande allowed region the bound by Chooz |J amounts to \U e s\ ^ 0.2. In the region 
not covered by Chooz |[/ e3 | can even be larger JTD|, |TT|. Thus neglecting \U^\ with respect to s 
in eq. (Q) is not really justified. Also note that in presence of a large hierarchy |m 3 | ^> |mi i2 | 
the effect of neglected parameters in eq. (U) can be enhanced by m 3 /mi j2 and produce seizable 
corrections. A non vanishing U e s term can lead to different {m v )\2 and (mju terms. Similarly 
a deviation from maximal mixing ^ U T 3 distorts the e 2 terms in the 23 sector of m v . 
Therefore, especially for a large hierarchy, there is more freedom in the small terms in order to 
construct a model that fits the data than it is apparent from eq. (|]). 

Given the observed frequencies and our notation in eq. there are three possible hierar- 
chies of mass eigenvalues: 



A 
B 

C 



|m 3 | >> |m 2 ,i| 

|mi| ~ |m 2 | >> |m 3 | 

|mi| ~ |m 2 | ~ |m 3 | 



(6) 



(in case A there is no prejudice on the mi, m 2 relation). For B and C different subcases are 
then generated according to the relative sign assignments for mi j2j3 . For each case we can set to 
zero the small masses and mixing angles and find the effective light neutrino matrices which are 
obtained both for double and single maximal mixing. Note that here we are working in the basis 
where the charged lepton masses are diagonal, and approximately given by m\ = Diag[0, 0, m T ]. 
For model building one has to arrange both the charged lepton and the neutrino mass matrices 
so that the neutrino results coincide with those given here after diagonalization of charged 
leptons. For example, in case A, mdiag = Diag[0, 0, m 3 ] and we obtain 



m u /m 3 










l 

2 







(7) 
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In this particular case the results are the same for double and single maximal mixing. Note 
that the signs correspond to the phase convention adopted in eq. (|2|). If one prefers all signs 
to be positive it is sufficient to invert the sign of the third row of the matrix U in eq. (H). We 
can similarly proceed in the other cases and we obtain the results in table I (where the overall 
mass scale was dropped), which we now discuss. 







double 
maximal 
mixing 


single 
maximal 
mixing 


A 


Diag[0,0,l] 


[0 1 

1/2 -1/2 
[o -1/2 1/2 J 




"0 

1/2 -1/2 
-1/2 1/2 




Bl 


Diag[l,-1,0] 




1/V2 1/V2' 
1/V2 
1/V2 






"10 

-1/2 -1/2 
-1/2 -1/2 




B2 


Diag[l,l,0] 


[10 1 

1/2 1/2 
|_0 1/2 1/2 J 




1 
1/2 1/2 
1/2 1/2 




CI 


Diag[- 1,1,1] 


[ -1/V2 -l/\/2] 

-1/V2 1/2 -1/2 
[-I/V2 -1/2 1/2 J 




"-1 0" 
1 
_ 1_ 




C2 


Diag[l,-l,l] 




1/V2 1/V2' 

1/V2 1/2 -1/2 
1/V2 -1/2 1/2 






1 
0-1 
0-10 




C3 


Diag[l,l,-1] 


[1 0] 
1 

|_o 1 oj 




"10 0" 
1 
1 





Table I : Zeroth order form of the neutrino mass matrix for double and single maximal 
mixing, according to the different possible hierarchies given in eq. (||). 

The completely degenerate case C is the only one that could in principle accommodate 
neutrinos as hot dark matter together with solar and atmospheric neutrino oscillations. For 
this the common mass should be around 1-3 eV. Then the solar frequency could be given by a 
small 1-2 splitting, while the atmospheric frequency could be given by a still small but much 
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larger 1,2-3 splitting. Note that we have not included in Table 1 the degenerate case C with 
three equal signs. Clearly the mixing matrix in this case is completely determined by the small 
degeneracy splitting corrections and the zeroth order approximation is irrelevant. A strong 
constraint arises in the degenerate case from neutrinoless double beta decay which requires 
that the ee entry of m v must obey |(mj,) ee | < 0.46 eV [T2|]. As observed in ref. |13|, this bound 



can only be satisfied if bimixing (that is double maximal mixing) is realized. In fact we see 
from the expression of e in eq. (]5|) that we need m x ~ — m 2 and c 2 ~ s 2 for a cancellation 
to occur in (m u ) ee . Since the solar mixing can only be either near maximal or very small, 
only the bimixing solutions CI and C2 survive (which are physically equivalent []), as can be 
verified in table 1. So all other solutions of type C can only survive if the common mass is 
below 0.46 eV. We think that it is not at all clear at the moment that a hot dark matter 
component is really needed ||14j| . However the main reason to consider the fully degenerate 
solution is that it is compatible with hot dark matter, so that only the solution C1/C2 in the 
case of bimixing is interesting. Note that for degenerate masses with m ~ 1 — 3 eV we need 
a relative splitting Am/m ~ Am 2 atm /2m 2 ~ 10~ 3 — 10~ 4 and a much smaller one for solar 
neutrinos explained by vacuum oscillations: Am/m ~ 10~ 10 — 10~ n . We are unable to imagine 
a natural mechanism compatible with unification and the see-saw mechanism leading to such 



a precise near symmetry. An ansatz of this sort has been proposed in ref. |15 



Of the remaining possibilities two are particularly remarkable for their simplicity: solution 
A ( e Q- (|ZD) an d solution Bl ( for bimixing). Of the six independent entries of a symmetric 
matrix three (four) are zero for solution A (Bl bimixing). However, solution A requires the 
entries in the 23 sector to be related (in particular the 23 subdeterminant must be zero, in order 
for two of the eigenvalues to be vanishing in zeroth order, as assumed in this case), solution Bl 
requires the two non zero entries in the first row or column to be nearly equal. In the following 
we shall discuss how these patterns can be generated in a natural way. 

As a first orientation we observe that, after diagonalization of the charged lepton Dirac 
mass matrix, we still have the freedom of a change of basis for the right-handed neutrino fields. 
In fact right-handed charged lepton and neutrino fields, as opposed to left-handed fields, are 
uncorrelated by SU(2) (g) U(l) gauge symmetry. We can use this freedom to make the Majorana 
matrix diagonal: M~ l = T duO with du = Diag[l/ri,l/r2,l/r^\. Then if we parametrize 
the matrix Omu by z^ we have: 

m vab = (m^M _1 m D ) a6 = V ZcaZcb . (8) 

c r c 

From this expression we see that, while we can always arrange the twelve parameters z a t> and 
r a to arbitrarily fix the six independent matrix elements of m u , case A is special in that it 
can be approximately reproduced in two particularly simple ways, without relying on precise 
cancellations among different terms: 

i) One of the right-handed neutrinos is particularly light and, in first approximation, it is only 
coupled to fi and r fllS] . Thus, r c ~ r\ (small) and z c \ ~ 0. In this case JD], |T7|] the 23 



subdeterminant vanishes, and one only needs the ratio \z c2 /z c3 \ to be close to 1. 
1 The solutions CI and C2 are related by the exchange s «-* c. 
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ii) There are only two large entries in the z matrix, \z c2 \ ~ | z c % | , and the three eigenvalues r a are 
of comparable magnitude (or, at least, with a less pronounced hierarchy than for the z matrix 
elements). Then, again, the subdeterminant 23 vanishes and one only needs the ratio | z^/z^ 1 
to be close to 1. 

The possibility ii), which was noticed in our previous paper ||, involves no more fine tuning 
than i). One solution is more specific on M, the other on mjj. Note that solution Bl requires 
a more delicate balance between the Dirac and Majorana matrices, which, however, can be 
realized in particular models [|ll|, |17j], as we shall see in the following. 



We now discuss the mechanism ii) in more detail. We start with the order zero approxima- 
tion, as in eq. (|7|) or table 1, where all small entries in the matrix are set to zero. Assume that, 
in the basis where charged leptons are diagonal, the Dirac matrix mp, defined by the bilinear 
^bT^d^l-, takes the approximate form: 



mo oc 






x 



(9) 



This matrix has the property that for a generic Majorana matrix M one finds: 



m T D M~ x 



mo oc 





x 2 

X 



(10) 



1 is that the 33 entry is non zero. The reason for this insensitivity to 
@ can be diagonalized by a transformation of the form V T mr,U = 



The only condition on M~ 
M is that ttid given by eq 
d oc Diag[0, 0,1] with V (the right-handed field transformation) given by a block-diagonal 
matrix with (V)^^ = 1. Then we can write m v = UdV T M~ l VdU T . The diagonal matrix d acts 
like a projector so that only the 33 entry of dV T M~ 1 Vd is non vanishing. Here it is crucial that 
the neutrino eigenvalues are approximately (0, 0,m3). Finally the U transformation (the left- 
handed field rotation) brings this matrix into the form of eq. (|Tt]) that for \x\ — > 1 approaches 
that required by solution A. Note that the 23 subdeterminant of m v is vanishing for all values 
of x. The matrix U is directly the neutrino mixing matrix. Its form is: 



with 





c 


—s 





u = 
















x/r , 




1/r 


r 



X A 



(12) 



We see that at \x\ 7^ 1, U e 3 still vanishes but the atmospheric neutrino mixing is no more 
maximal. Rather one has 

sin 2 29 = 4sl£= „ AX \.„ . (13) 



77 



(1 + x 2 ) 2 

Thus the bound sin 2 29 ^ 0.8 translates into 0.6 ^ |x| ^ 1.6. It is interesting to recall that in 
ref. |Tj| it was shown that the mixing angle can be amplified by the running from a large mass 
scale down to low energy. 
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A zeroth order texture as in eq. @ with generic entries (m^) 2 3 = A and (7710)33 — B 
can well apply for both charged leptons and neutrinos. If for charged leptons we start with 
Ai and Bi and for neutrinos with A v and B u , after diagonalisation of the charged leptons by 
a left-handed matrix of the form in eq. ( |Tl~D with s — 0, one obtains a matrix of the same 
form with A/B = x, as required for mj) in eq. Precisely the value of \x\ is given by 

|x| = \{AiB u - BiA v )/ (A L A U + BiB v )\. 

It remains for us to discuss if the proposed strategy is not completely ruined by going beyond 
the zeroth order approximation. In fact terms of order A, A 2 ... have to be added to mn to 
introduce a hierarchical structure as observed for quarks and leptons. Since, in general, also M 
will be hierarchical in terms of some other small parameter rj, M _1 can contain terms of order 

I/77, 1/rf In the see-saw product, the ratios of a power of A and a power of 77 can generate 

unwanted additional terms of order one or even larger that can spoil the zeroth order result. 
Obviously a necessary condition is that the Majorana hierarchy is not too strong in comparison 
of the Dirac hierarchy, or that rj is not smaller than some power of A. We discuss this issue by 
showing that reasonable models exist where the mechanism is indeed realised. 

An attractive approach to reproduce the observed features of quark and lepton mass matrices 
is in terms of horizontal symmetries, in particular, for simplicity, U(l) charges ffjTsfl . When the 
charges do not match by n the corresponding mass term is suppressed by n powers of a small 
parameter that arises from the vacuum expectation value (divided by some large mass) of one 
or more extra scalar fields 9j that carry one unit of charge. There can be two different fields 
with opposite charges 6± f| so that the suppression factors can involve two different small 
parameters, according to the sign of n. We now show that our mechanism can be realised in a 
simple U(l) model with two oppositely charged fields 6±. 

We want to construct an example of charge assignments that lead to the desired form for the 
Dirac and Majorana matrices. In particular, we demand that, for the Majorana matrix M, the 
determinant is of order 1 (because we want to avoid that one lightest eigenvalue exists, which 
would be typical of mechanism i), while our case prefers comparable eigenvalues) and the same 
is true for the matrix element (M -1 ^ (which plays a special role in our mechanism). One 
possibility is to start with the following charges for the lepton doublet L and the right-handed 
neutrino fields: 

L = (2,0,0); it = (1,-1,0) (14) 

which lead to the textures: 





"A 3 


A 


A" 




"A 2 


1 


A 


mo oc 


A 


A' 


A' 


; M oc 


1 


A' 2 


A' 




A 2 


1 


1 




A 


A' 


1 



(15) 



It is simple to verify that, for generic coefficients, the effective light neutrino matrix is of the 
form: 

"A 4 A 2 A 2 ' 

m„ x A 2 1 1 , (16) 
A 2 1 1 



2 We regard the conjugate field 9* as distinct from 6. Therefore 9_ may coincide with 
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However, the crucial property is that the 23 subdeterminant of m u automatically vanishes in the 
limit A — ► 0, because of the simple mathematical property displayed by eqs. (|), ([3~0|) . Actually, 
for generic coefficients, the 23 minor is of order AA'. We now specify the charges of right-handed 
charged leptons. The requirements are that the hierarchy of eigenvalues corresponds to the 
observed e, /i, r mass ratios and that after diagonalisation of the charged leptons the form of 
the Dirac matrix for neutrinos is not spoiled (for this we need that the left-handed diagonalising 
matrix is of the form in eq. (|TTD with s small). An acceptable choice for the right-handed lepton 
charges is: 

^ = (-4,1,0) (17) 



which leads to: 

mi oc 



A' 2 A' 4 A' 4 
A 3 A A 
A 2 1 1 



This matrix has eigenvalues of order (A' 2 , A, 1). For A ~ A' ~ 1CT 1 — 1CT 2 and generic coefficients 
one obtains numerically reasonable solutions. In this analysis we have also taken into account 
the most general kinetic terms compatible with the chosen charge assignement. We found that 
numerically acceptable solutions are generated even in the presence of these terms. (See also 
ref. lyj). Most of them are of the bimixing type. This is because, as seen from eq. ([T6|), the 
11 entry of m u is much smaller than the 12, 13 entries. Given the general expression eq. (^) 
this particular pattern corresponds to double nearly maximal mixing. 

It is interesting to note that if all the U(l) charges for the left-handed lepton doublets and the 
right-handed anti-neutrinos have the same sign, the relation between the Dirac and Majorana 
hierarchies is too tight for the mechanism to work. In fact in this case, the neutrino Dirac and 
Majorana mass matrices can be written as |30| m v D = Q u rY d Q l and M = Q u rYmQvR where 
the Q's are diagonal matrices of the form Q = Diag(X qi , A 92 , A 93 ) in terms of one single small 
parameter A and the charges g^. Y are generic non hierarchical Yukawa matrices. It is easy to 
see that in constructing m u the right-handed charges drop away and m v = Q lY^Y^ 1 Y d Q L . 
With the above left-handed charge assignments the resulting matrix m u indeed has the 22, 
23, 32, 33 entries of order 1 while all others are suppressed, but the 23 subdeterminant is not 
vanishing because the matrix Y^Y^Yr, is completely generic. Thus the resulting eigenvalue 
structure is not the right one in general and must be fixed by hand. It is easy to see why our 
mechanism does not work in this case. In fact, the same reason that ensures the cancellation 
of right-handed charges in the see-saw product makes the connection between the hierarchies 
of the Dirac and Majorana matrices such that the A terms of mo and the 1/A terms of M _1 
necessarily generate finite terms that spoil the vanishing of the 23 subdeterminant. So models 
of type A cannot in general be explained in terms of all positive (or all negative) L and R v 
charges. For our mechanism to work we need that the Dirac and Majorana mass matrices are 
generated by less tightly connected terms, as in the above example. This could also be the 



case in models based on GUT's p2| , for example Susy SO(10) or flipped SU(5), where mn 
and M are generated by different Higgs multiplets and representations. Note that the model 
Bl (bimixing) can be reproduced in terms of left-handed charges (a,-a,-a) for (eL, ^l,tl), as 
observed in ref. [O], |T7|, but the 12 and 23 entries, of order 1, must be arranged to be nearly 
equal. 
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In conclusion, we have studied a mechanism, introduced in |J , that can explain the observed 
pattern of neutrino oscillations, both for double and single maximal mixing. In our case both for 
Dirac and light neutrino masses there is one dominant eigenvalue: |m3| >> 1 771-1,2 1 as m case A 
of table 1. Our proposed mechanism only contains a minimum of fine tuning and compares well 
with all other mechanisms sofar presented |22|, ^3|, [24], ^5], |26 



Some of them were also discussed 
here. From these examples we learn that after all large mixings are not so unplausible as it 
appears at first sight. In all of these models there is no reason why the atmospheric neutrino 
mixing should be maximal but it can be large. In many cases once the atmospheric mixing is 
large also the solar mixing turns out to be large, but this is not necessarily true in all models, 
for example in our model. 
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